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Abstract. We study the role of impact parameter on the collective flow and its disappearance 
for different mass asymmetric reactions. The mass asymmetry is varied from to 0.7 keeping 
the total mass of the system fixed. Our results clearly indicate a significant role of impact 
parameter on the collective flow and its disappearance for the mass asymmetric reactions. The 
impact parameter dependence is also found to vary with mass asymmetry of the reaction. 



1. Introduction 

On the collision of two nuclei, the density and pressure increases in the interaction region and 
at finite colliding geometries, there is an inherent asymmetry in the pressure which results in 
a collective transverse flow of matter toward the direction of lowest pressure. This collective 
transverse flow is one of the most extensively used observable to study the equation of state (EoS) 
as well as in- medium nucleon-nucleon (nn) cross-section of the nuclear matter. At the particular 
value of energy where the attractive scattering (dominant at energies around 10 MeV/nucleon) 
balances the repulsive interactions (dominant at energies around 400 MeV/nucleon), collective 
transverse flow in the reaction plane disappears. This energy is termed as Balance Energy (E(, a i) 

[H El El HI El El El El ED]. 

Both collective flow and E^ai are found to be highly sensitive towards the nn cross-section 
[SSHOElElElEldO], size of the system (A TO t = A T + A P ; where A T and A P are the masses 
of the target and projectile, respectively) [4j, mass asymmetry of the reaction (77 = ^ T ~^ p ) [TT] . 
nuclear matter equation of state [21 El [3J El El El El [TO] , incident energy (E MeV/nucleon) [12] 
as well as the colliding geometry (b = j^—; where b max = Ri + R2; Ri is the radius of projectile 
or target) [21 El El El UD]- Several experimental and theoretical attempts have been employed 
in order to explain and understand these observations [IJE1E1I2E1ISC3E1E1ITO1III1CE2]- Due 
to the decrease in the compression reached in heavy-ion collisions with increase in the impact 
parameter, E^ a i is found to increase approximately linearly as a function of the impact parameter 
|13t fl"4"] . But most of the studies in the literature mainly focus on the symmetric and nearly 
symmetric reactions. Our present aim, therefore, is at least twofold. (1) To study the variation 
of collective transverse flow with impact parameter for different r\ keeping the incident energy 
fixed. (2) To study the sensitivity of collective transverse flow and its disappearance on the mass 
asymmetry of the reaction at different impact parameters. The Quantum Molecular Dynamics 
(QMD) model [S1[TOJ[]I>] used for the present analysis is explained in the section 2. Results and 
discussion are presented in section 3 followed by summary in section 4. 



2. Description of the model 

In the QMD model, each nucleon propagates under the influence of mutual two- and three-body 
interactions. The propagation is governed by the classical equations of motion: 



where the Hamiltonian is given by 



dt dpi ' 

dpi _dH_ , , 

dt dn' U 
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Our total interaction potential V tot reads as [HJ [151 US] 

ytot yLoc _|_ -yYuk _j_ -yCoul 

with 

= /^(f- - rj) + t 2 5(f- - f,,)^ - r fc ), (5) 
= tge-lfi-^l/m/ ( | rl _ fj | /m) ; (6) 

with m = 1.5 fm and t$ = -6.66 MeV. 

The static (local) Skyrme interaction can further be parametrized as: 

Here a, (3 and 7 are the parameters that define equation of state. The parameters a, /3, and 
7 in above Eq. (7) must be adjusted so as to reproduce the ground state properties of the 
nuclear matter. The set of parameters corresponding to different equations of state can be 
found in Ref. [IB]- It is worth mentioning that as shown by Puri and coworkers, Skyrme 
forces are very successful in the analysis of low energy phenomena such as fusion, fission and 
cluster-radioactivity, where nuclear potential plays an important role |17j . 

3. Results and Discussion 

For the present study, we simulated the reactions of f^Ne +10 Ne (77 = 0), g 7 +ff Na (77 = 0.1), 
fN+HMg (77 = 0.3), fB+USi (77 = 0.5), and ^Li+^S (77 = 0.7) for total mass (A TOT ) = 40, 
f Ca +|° Ca (77 = 0), f 8 Ar +% Ca (77 = 0.1), $Si +f A Cr (77 = 0.3), j° Ne +«j Ni (77 = 0.5), and 
fB + 7 3 2 Ge (77 = 0.7) for total mass (A TOT ) = 80, f 6 Kr +fg Kr (77 = 0), l%Ge +|g Zr (77 = 0.1), 
liFe +1^ Cd (77 = 0.3), \lCa +g° Te (77 = 0.5), and \\Mg +lf Ce (77 = 0.7) for total mass 
(A TOT ) = 160, and ™>Te +\f Te (77 = 0), \™Cd +lf Ba (77 = 0.1), gjSr +66 6 Dy (77 = 0.3), 
llNi +\f W (77 = 0.5), and f^Ar Pb (77 = 0.7) for total mass (A TOT ) = 240. The impact 
parameter is varied from b/b max = to 1 in small steps of 0.25. The charges are chosen in a 
way so that colliding nuclei are stable nuclides. A soft equation of state with isotropic energy 
dependent cugnon cross-section (labeled as Soft* 50 ) is used for the present calculations. 

The balance energy (E), a i) is calculated using the directed transverse momentum < P x ir >, 
which is defined as: 

{Pf r ) = jT, s[ ^{Y({)}p x (i), (8) 
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Figure 1. (Color online) The < P x ir > (MeV/c) as a function of reduced impact parameter 
(b/b max ) for different system masses. The results for different mass asymmetries ry = 0, 0.1, 
0.3, 0.5, and 0.7 are represented, respectively, by the solid squares, circles, triangles, inverted 
triangles, and diamonds. Results are at an incident energy of 200 MeV/nucleon. 



where Y(i) and p x (i) are the rapidity distribution and transverse momentum of i th particle, 
respectively. 

In Fig. 1, we display at a fixed energy, the < P x ir > as a function of reduced impact 
parameter (b/b max ) for rj = 0-0.7, keeping total system mass fixed as 40, 80, 160, and 240. 
All reactions are followed till 200 fm/c, where < P x ivr > saturates. In all cases < P x %r > first 
increases with increase in impact parameter, reaches a maximal value and after passing through 
a zero at some intermediate value of impact parameter, attains negative values. The trend 
is uniform throughout the mass asymmetry range. The value of impact parameter at which 
< P x %r > attains a zero (which is termed as Geometry of Vanishing Flow (GVF) [18]) varies 
with r\ and Atot- For lighter systems and larger r/, the value of GVF is smaller compared to 
the heavier systems and smaller rj. 

In Fig. 2, the variation of GVF as a function of rj is displayed for different system masses. 
The percentage variation in GVF while going from 77 = to 0.7 is -40.9%, -24.59%, -16.44%, 
and -20.93%, respectively for Atot = 40, 80, 160, and 240. It is clear from the figure that the 
effect of mass asymmetry of the reaction on GVF decreases with increase in system mass. This 




Figure 2. (Color online) The geometry of vanishing flow (GVF) as a function of 77 for different 
system masses. The results for different system masses (Atot) = 40, 80, 160, and 240 are 
represented, respectively, by the half filled squares, circles, triangles, and inverted triangles. 



is similar to as predicted for E^i [TT] . Due to the decrease in nn collisions with increase in 77 and 
impact parameter and increase in Coulomb repulsion with increase in Atot, the Ef, a i increases 
with increase in 77 and impact parameter, while it decreases with increase in Atot- Since the 
present study is at fixed incident energy, therefore, the value of impact parameter, where flow 
vanishes, decreases as 77 increases. 

In Fig. 3, we display < P^ ir > as a function of incident energy ranging between 40 
MeV/nucleon and 800 MeV/nucleon, for different mass asymmetric reactions keeping the total 
mass of the system fixed as 240. The results are shown for different impact parameters. In 
all the cases, transverse momentum is negative at lower incident energies which turns positive 
at relatively higher incident energies. The value of the abscissa at zero value of < P^ tr > 
corresponds to the balance energy. The figure indicates that (i) for all values of 77 and impact 
parameters, the transverse momentum increases monotonically with increase in the incident 
energy. The increase is sharp at smaller incident energies compared to higher incident energies 
where it starts saturating, (ii) due to decrease in overlap volume and hence number of collisions 
with increase in 77, the transverse momentum starts suppressing as 7/ increases and hence the 
balance energy increases with increase in r\. (iii) the variation in < P^ ir > with 77 decreases 
with increase in impact parameter. But with increase in impact parameter the balance energy 
increases for all values of 77. 

The above mentioned findings can be understood by decomposing the total transverse 
momentum into contributions due to mean field and two- body nn collisions as shown in Fig. 4. 
The symbols are explained in the caption of the figure. One notices that mean field flow increases 
up to couple of hundred MeV/nucleon and then saturates. In the lower incident energy region, 
the flow due to mean field is smaller for larger asymmetries compared to smaller asymmetries. At 
higher incident energies the mean field flow is nearly independent of 77. Similar behavior is seen 




Figure 3. (Color online) The < P^ ir > (MeV/c) as a function of incident energy for system 
mass Atot = 240. The results are shown for different mass asymmetries (77 = 0-0.7) at reduced 
impact parameters b/b max = 0.25 (upper panel), 0.5 (middle panel), and 0.75 (bottom panel). 
The lines are only to guide the eye. Symbols have the same meaning as in Fig. 1. 
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Figure 4. (Color online) The decomposition of < P~ ir > displayed in Fig. 3 into mean field 
(half filled symbols) and collision part (open symbols) as a function of incident energy. 
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Figure 5. (Color online) The £7& a ; as a function of r/ for system mass Atot = 240. The results 
for different impact parameters b/b max = 0.25, 0.5, and 0.75 are represented, respectively, by 
the crossed squares, circles, and triangles. Lines are to guide the eye. 



for other impact parameters also. The flow due to binary nn collisions, however, increases with 
increase in incident energy for all 77, but as the impact parameter increases, the nn collision flow 
starts saturating. Obviously as explained earlier, the flow due to binary nn collisions decreases 
with increase in 77. 

In Fig. 5, we display E^ai as a function of tj for b/b max = 0.25, 0.5, and 0.75, keeping Atot 
fixed as 240. Various symbols are explained in the caption of the figure. From the figure, we 
see that Eb a i increases with increase in tj for all values of b/b max - Also the 77 dependence of E^i 
increases with increase in impact parameter. 

In Fig. 6, we display the percentage change in balance energy AE ba b r ax (%), defined 

as A£^ 6mo *(%) = ((E^ ma ' ¥0 - 2S -E^ ma ' =0 ' 25 )/EQ mam=0 ' 25 )xlOO as a function of 77. Lines 
represents the mean value of variation. It is clear from the figure that the effect of impact 
parameter variation is almost uniform throughout the asymmetry range for every fixed system 
mass. We also found that the mean variation decreases with increase in Atot- Clearly, the 
effect of impact parameter variation is independent of 77. 

In Fig. 7, we display the percentage difference AE^ al (%) defined as AE^ al (%)) = {(E^ - 
Ebai°)/Ebai°)>< 100 as a function of reduced impact parameter (b/b max ). Lines are the linear 
fits (oc m , b ). It is clear from the figure that the effect of the asymmetry variation increases 

Umax 

with increase in the impact parameter for each mass range. This is due to the fact that with 
increase in impact parameter, the number of binary nn collisions decreases and the increase of 
mass asymmetry further adds the same effect. 
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Figure 6. (Color online) The percentage difference AE^ max (%) as a function of r\ for different 
system masses. The results of the percentage difference for different colliding geometries b/b max 
= 0.5 and 0.75 are represented, respectively, by the crossed circles and triangles. Horizontal 
lines represent the mean value of AE^J max (%) for each b/b max . 



4. Summary 

We presented a detailed study on the role of impact parameter on the collective flow and its 
disappearance for different mass asymmetric reactions using the quantum molecular dynamics 
model. For the present study, the mass asymmetry of the reaction is varied from to 0.7 
by keeping the total mass of the system fixed. A significant role of impact parameter on the 
collective flow and its disappearance for the mass asymmetric reactions is seen. The impact 
parameter dependence is also found to vary with mass asymmetry of the reaction. 
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Figure 7. (Color online) The percentage difference AE^ al (%) as a function of b/b max for different 
system masses. The results of the percentage difference for different asymmetries r/ = 0.1, 0.3, 
0.5, and 0.7 are represented, respectively, by the solid circles, triangles, inverted triangles, and 
diamonds. Lines are the linear fits (oc m b b ); m values without errors are displayed. 
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